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Abstract. We consider the spectrum of two ultracold harmonically trapped atoms 
interacting via short-range interactions. The Green's function approach is used to 
unify the two and three dimensional cases. We derive criteria for the universality of 
the spectrum, i.e. its independence of the details of the short-range interaction. The 
results in three dimensions are examplified for narrow s-wave Feshbach resonances and 
we show how effective range corrections can modiiy the rearrangement of the level 
structure. However, this requires extremely narrow resonances or very tight traps that 
are not currently experimentally available. In the two-dimensional case we discuss the 
p-wave channel in detail and demonstrate how the non-universality of the spectrum 
arises within the Green's function approach. We then show that the spectrum is 
not particularly sensitive to the short-distance details in the case when the two-body 
interaction has a bound state. 
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1. Introduction 

In an age of rapidly increasing computational power, exact methods and benchmark 
solutions continue to have tremendous importance as a means of gauging numerical 
calculations and provide invaluable analytical insights [1]. Ultracold atomic gases have 
emerged as a field with great potential as a laboratory benchmark for many fields of 
physics. The extreme control exercised over the systems in terms of trapping geometry 
and inter-atomic interactions allow experimenters to prepare samples that simulate the 
intricacies of many different models that are applied in other fields of physics and the 
term 'quantum simulator' is often applied [2\. The clean conditions provide hope of 
understanding some of the paradigmatic models such as the Hubbard model that are 
used frequently in the study of materials and other condensed-matter systems. 

In the case of cold atoms, the interactions are typically short-ranged and the samples 
have to be in an external confinement, most often provided by a magnetic or optical 
potential. In the case of a harmonic oscillator trapping potential, it turns out that 
the problem of two atoms interacting through a short-range potential can be exactly 
solved as demonstrated by Busch et al. [3]. The authors of Ref. [3] approximated the 
two-body interaction potential by a zero-range pseudopotential introduced long ago by 
Hellmann [5] and Fermi [5] . The predictions of this exact model was subsequently tested 
experimentally in an optical lattice and found to be a very accurate description of the 
two-atom system [6j. 



The pseudopotential approach has become somewhat of a paradigm itself in cold 
atoms. Its success can largely be attributed to a seperation of scales; the two-body 
collisions energy is small at the low temperatures one usually aims for, and the density 
of the system is also much lower than typical matter densities. One can then model 
the interaction using only a few low-energy parameters such as the scattering length, 
a, and effective range, rg [7]. The true range of the potential, given by the van der 
Waals length for neutral atoms, is much smaller than the interparticle spacing. The 
parameters a and r^ therefore characterize the system, and since these are independent 
of the shape of the two-body potential, one refers to this as a universal regime, i.e. 
when |a|,n~^/^ ^ rg with n the density. A really nice feature of the cold atomic gas 
system is the tunability of the interaction parameters through Feshbach resonances [8] , 
which allows one to explore the full region of parameter space, including the interesting 
unitarity limit where l/|a| — )■ 0. 

In the case of a trapped system with harmonic oscillator length, 6, a fundamental 
question concerns the binding energy of an A^-body system in the universal regime 
when we also require that ro <C b. For equal mass particles, this has been studied 
numerically using a host of different methods [9]. Interestingly, for a — oo, the three- 
body problem in a trap can be exactly solved [lO]. In the same spirit, exactly solvable 
models in a harmonic approximation approach have been proposed [11]. Very recently, 
experiments in so-called microtraps have demonstrated that few-body systems can in 
fact be produced with cold atoms and universality and shell structure can be explored 
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[12] • The general framework of effective field theory is very suitable for problems with 
seperation of scales, and it has been applied successfully to the three-body problem in 
cold atoms [13]. Recently, there has been a lot of interest in applying these techniques 
within an oscillator basis [HI [15] to address few-body Fermi systems in cold atoms 
and in nuclei within the no-core shell-model approach [16]. Similar methods have also 
been used to study few-body bosonic systems in traps (see [T7| and [18] for details and 
references) . 

The model of Busch et al. is the basic foundation upon which many of the 
developments discussed above reside. Here we take a fresh look at the model from 
a Green's function point of view [211 l20l 119] . This is done within a two-channel 
formalism. We consider the general angular momentum /-wave case and derive the 
expression for the eigenspectrum in three dimensions as a function of the scattering 
phase shift. For I > 2, there are obstructing terms that imply a dependence on the 
short-range two-body potential, i.e. the spectrum is non-universal. A criterion for 
the applicability of the universal formula is subsequently derived. As an example, we 
consider a two-channel model for narrow Feshbach resonances and show that interesting 
spectral changes occur when including the effective range term. This can be interpreted 
in terms of the Zeldovich rearrangement effect [22] which occurs in systems with a long- 
range attractive potential and a short-range two-body attraction that dominates at 
small distance (some recent discussion of the effect can be found in Refs. [231 El])- The 
observation of these effects in experiments require, however, extremely tight trapping 
potentials or extremely narrow Feshbach resonances, both of which are beyond the 
current experimental capabilities. 

Two-dimensional setups are currently of great interest in the cold atom community 
and a number of impressive experimental results have been reported recently [25j . 
In the second part of the paper we consider a two-dimensional geometry and derive 
the eigenspectrum within the Green's function approach. We find that for angular 
momentum m > 1, there are non-universal terms, i.e. a spectrum that depends on the 
short-range details of the two-body potential. Emphasis is put on the p-wave case m = 1, 
where we find a closed formula for the spectrum in terms of the low-energy parameters of 
the interaction. To obtain the eigenspectrum, we use a generic form of the p-wave phase- 
shift in two dimensions, which is similar to that obtained in hard-sphere or square well 
potential models. The p-wave spectrum in the case where the interaction allows a two- 
body bound state is very similar to that of s-waves. The shape of the p-wave spectrum 
is almost universal, depending only slightly on the exact potential model used. 

2. Basic Two-Channel Formalism 

The physics around Feshbach resonances [8J is most naturally by models which 
explicitly take the open (scattering) channel and the closed (molecular) channel into 
account [26l |8] . Here we consider such a model within the Green's function approach 
[2T] . The setup has a dressed state with open, \E'(r), and closed channel, /30(r) 
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wavefunctions, where /3 is the amphtude of the closed channel and the normalization is 
+ J (ir|\l/(r)p = 1 since we assume that is normalized. The Hilbert space of the 
closed channel is therefore one- dimensional for simplicity. Notice that (f){r) has angular 
momentum / and projection m/. The wave eqautions are 

Vtrap^ir) + W{r)/3<j){r) = E^{r) (1) 
W{r)^{r) + E^Mir) = EI3(t>{r), (2) 

where the operator is Vtrap = T)free + Vix)- Here V{r) is the trapping potential which 
we assume to be an isotropic harmonic oscillator with trap length h = ^Jfl/ where 
u is the oscillator frequency and fi is the reduced mass, i.e. V{r) = ^^u'^r'^ The free 
particle operator has the standard form Vfree = —h^\/2ii. The energy of the closed 
channel molecule is denoted Emoi and the coupling between the channels is parametrized 
by the real function W{r). We assume that the range of W{r) is much smaller than h. 



3. Three Dimensions 



The 3D Green's function GE{r,r') is defined by 



If we define 



[V,rap - E] GE{r, r') = ^6{r - r'). (3) 

A* 

:F = E-E.^,i + -^I drdr'^*{r)GEir, r')0(r'), (4) 

where 0(r) = W{r)(f){r), then upon substitution of Eq. [3] into Eqs. [1] and [2] we obtain 
= 0. In momentum space we have 0(fc) = \fATii^k^Yim{k)aim-, which defines the 
coupling constant 



aim = V4^k-' j dr4){r)]i{kr)Yi*^{r). (5) 

In order to relate the spectrum in the trap to the free-particle scattering properties 
we have to calculate the scattering amplitude. Therefore we consider the scattering 
problem 

"^{r) = e^'^- - ^(3 I dr'4>ir')G%ir, r'), (6) 
where G^{r,r') is the free-particle Green's function and E = h?k'^/2fi is the scattering 



energy (see [Appendix A|. Solving for /3 in Eq. [2]and inserting into the scattering solution 



in Eq. [6l we find 



^' r(k)fdr'4,{r')G%(r,r') 



where J-q is defined analogous to J-" above but with the free Green's function G^{r, r'). 
The solution of Eq. [3] for out-going wave boundary conditions when the trap is absent 
has the asymptotic behavior 

„ik\r—r'\ „ikr 

G'0(r,r') = r e"'*''"'' for r ^ cx), (8) 



Two-Body Spectra of Harmonically Trapped Atoms 

where k' is the final momentum which fulfills |fc| = Using the fact that 

we finally obtain 

fk{k') : _ 



Through partial wave decomposition fk{k') = 4:7iJ2im fi.m{k)Yi^{k)Yim{k') we get 



dre-'''-''(j)*{r), 
^~^*m{k') 



(9) 



(10) 



flm{k) 



(Emol — E) 

drdr'4)*{r)G%{r,r')4){r'). 



Since both G%{r, r') and G^{r, r') are singular at r = r', we have to regularize by 
isolating the finite part through (r, r') = G^{r, r') — G%{r, r'). We find 



\air. 



flmik) 



(12) 



(13) 



drdr'(f)*{r)G§{r,r')(f){r'), 

where G^{r,r') satisfies 

(Vtrap - E) Gi(r, r') = -V^(r)G°;(r, r'). 

By symmetry we only need to consider r > r', and we therefore write 

Gliry) = giir,r%ikr')Yilir)Y,Jr'), r > r'. (14) 

The solution for gi{r,r') can easily obtained by noting that a particular solution to 
Eq. [T3]is —G%{r,r') [19]. Adding the homogenoues solution gives 



gi[r,r 



A(r')(-)'M(-/ii,/ + 3/2;rV&') 





+B{r'){-y+'M{-f,2A/2-l;ryb') 
r 



-Arrikh'-pikr), 



(15) 



with Hi = ^ — 1/2 — 3/4 and /i2 = + 1/2 — 1/4. Here M(a, b; x) is the confiuent 
hypergeometric function and h\^\x) is the spherical Hankel function of the first kind. 
Demanding that G^{r, r') vanish as r — i- oo yields 

r(i/2-/)r(-/xi) 



A{r') 



-B{r'). 



(16) 



r(/ + 3/2)r(-/i2) 

Furthermore, by demanding that G^{r,r') be regular as r, r' — t- (while maintaining 
the condition r > r'), we demand that 

r(/ + i/2)2'+2y^ 



lim B{r') - ni.\i ■ 

r'^o b {kby 

gi{r, r') can now be determined for small r. For r — )■ and / < 2 

47r(A;r)' 



g[r,r 



(2/ + 1)!!A;2' 



{2l + iy.\A{r')k^ _ 
AtxU 



(17) 



(18) 
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where the leading term ji{kr) for A;r ^ 1 has been isolated to resemble the structure 
of Eq. [5l The imaginary part of Eq. [TSl comes directly from Im h[^\kr) . For I > 2 
there are additional terms at order with a < I. Inserting the solution into Eq. [12] and 
assuming that is short-ranged, we find 



(2/ + l)!!A(0)fc' 
A-kU 



ik 



21+1 



(19) 



where aim cancels which was the object of the regularization procedure. Notice that 
this does not depend on m due to the isotropy of the trap potential. Using Eqs. fT6land 
[T7t we finally have 



flrajk) 



\/+lo2/+l 



r( 



-1 



(20) 



r(-/i2) 

This has to be related to the scattering amplitude in terms of the /-wave phase-shift 
Si{k), which is 



flmik) 



P'+i cot 5(A;) - zP'+i 
Matching with the solution above we obtain 



_E_-' 



U 2 2hi^' 



22ZH 



{khf+^ cot 5{k), 



(21) 



(22) 



which recovers previous results [21 HSl Furthermore, the expression in Eq. [2^ shows 
how to include higher order terms from effective-range expansions. 

The result in Eq. [19] holds for short-ranged (p and for I < 2. There are generally 
terms obstructing the simple formula derived above for / > 2. These have order from 
a = 3 — /toa = / — lin steps of two up to the leading A{0) term of order /. For 
instance the / = 2 case has a term proportional to r, the / = 3 case has a constant and 
an term and so forth. This can be seen by considering the series expansion of h[^\x). 
Since we demand that (p is very short-ranged, we assume r/6 ^ 1. In general the most 
important term near r = will be r'^~' for I > 2, thus it diverges for / > 3 which must 
be compensated by the behavior of 0. A sensible criterion for this to happen is that 
A{0) should dominate over the most divergent obstructing term. For general / > 2, this 
gives us the inequality 

J oi /^rti , 1/2) /rN3- 



, ,r\ ,2'v^r(/ 



(23) 



(2/-3)(A;6)' 

W{r) is assumed to be of short-range and 0(r) is an /-wave wave function. Thus there 
will be some intermediate region around r = tq where has its weight. We want the 
condition to be satisfied at this distance. Inserting A{0) and using Eq. [22] we arrive at 

2^-3 21 + 1 [{21 -ly.lfl tan 6i{k)\ 



bJ 



> 



(24) 



2/ -3 4 (A;6)2'+i ' 

which agrees with Ref. [2^ for / = 2. The distance ro is essentially the same as the 
matching distance of Ref. [28] and Ref. [27] . We still require vo/b -C 1, so when the 
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(a) 



(b) 



-10 -5 



(c) 



(d) 



Figure 1. (Color online) 3D s-wave spectrum with afg/b = 0.01 and (a) r^g/b — 0.1, 
(b) r|o/6 = 1, (c) r|o/6 = 10, and (d) rfo/& = 100. The solid (black) lines are the full 
solutions while the dashed (red) lines are for r|g = 0. The dot-dashed (blue) horizontal 
lines indicate the non-interacting level structure (visible in the top left corner of (a) 
and (b)), while the dotted (blue) vertical lines are asymptotes Ihui above (or below) 
the non-interacting levels. 



right-hand side blows up, the universal formula breaks down and the details of the 
two-body interaction become important. 

For concreteness, let us consider a model for the coupling between the channels 
where W{r) = Woe~^^"''^ with oq the Bohr radius [29]. Assuming that 0(r) oc r', the 
maximum of (p is at tq = laQ. The left-hand side of Eq. [2l] is therefore very small 
as 6 ~ 1/im for typical traps. The factors depending on / on the right-hand side are 
increasing but only gradually. Using the lowest-order in the effective-range expansion, 
the criterion is (a//6)^'+^ ^ 1 where ai is the /-wave scattering length (whenever it is 
well-defined). Thus we conclude that the formula works only away from resonance. 

3.1. 3D Feshbach Model 

The formula in Eq. |22] has been discussed by a number of authors for both s-waves 
[31 [301 EH], p- waves [201 IlH], and waves [27]. In the context of Feshbach resonances, 
most studies have assumed that the effective range can be neglected which is true for 
wide resoanances [H]- To complement this, we considered here the opposite limit of 
very narrow resonances where the effective-range is large. This is interesting in cold 
atomic gases with two or more different species of atoms which typically have narrow 
resonances. In particular, at the point where the scattering length goes to zero we 
expect the corrections from higher-order terms to become important [31j. The universal 
behavior of one-channel models is described through one parameter, the scattering 
length. As a function of applied external magnetic field we parametrize the /-wave 
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field dependnent scattering length ai{B) in the following way 

AB 

''bg 



1 - 



(25) 



B-Bo_ 

where a\g is the background scattering length away from the resonance located ai B = Bq 
with width AB. Since we are interested in going beyond the one parameter description, 
we use a two-channel model of a Feshbach resonance [32] with corresponding open-open 
channel s-wave T-matrix 

= 1 , (26) 

where A/x is the difference in magnetic moment of the open and closed channels. 
Combining this with Eq. |20l we get 

Introducing the background effective range rfg = —ti?/{Ajj,AB^alg), and the useful 
quantities x = {B - Bo)/AB and f{E) = v^r(| - ^)/T{\ - we have 

1+ 2b^ \ I =l^/(^)- (28) 

We note that since a[gAfiAB > for all resonances with any / [8], we have r^Q < 
always. Isolating x yields 

2 b b hoo i_^f(^E)' 

In the case of a wide resonance or a large trap, Irfgl/fe — )■ 0, we recover ^^i^ = /(-E')- 

There is one caveat that has to be addressed before we proceed to study the 
effective range corrections to the two-body spectrum. This is related to the limit when 
the scattering length becomes where small. Here it is not necessarily clear that the 
properties of the spectrum will be universal in the sense that higher-order parameters 
from the effective range expansion (beyond a and re) can be neglected. This issue has 
been discussed in the context of effective field theory in Ref. |33j . 

The details of Feshbach resonances when the scattering length goes to zero has 
been considered for both trapped bosons and fermions in Ref. [31]. There it was found 
that the effective interaction is quadratic in the relative momentum at lowest order, 
since the usual constant piece proportional to a"^ vanishes. The coefficient depends on 
the background parameters of the resonance through the combination (ajg)^reo. The 
studies in Ref. [31] demonstrated that no anomalous behavior is seen when approaching 
the zero-crossing of the Feshbach resonance. 

Of course, if this quantity happens to be very small, even higher-order terms in 
the effective-range expansion must be taken into account. Here we are assuming that 
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Figure 2. (Color online) Same as Fig. [T] for al^/b = 0.1 and (a) r^p/fo = 0.1, (b) 
r|o/6 = 1, (c) r|o/6 = 10, and (d) = 100. 



kfegl ^ '^vdW niuch larger than the true range of the potential, given by the van 
der Waals length, ''"vdW' Likewise, for the narrow resonances we are interested in here 
I'^eol ^ '^vdW' ^^^^ regime that we expect the behavior to be universal, since 

this implies that low-energy scattering still dominates the two-body collisional dynamics. 
The smallest values used below are alg/b = 0.01. For a typical trap with b ~ l/xm, this 
is larger than T'vdW most atoms used in cold gas experiments. 

In Figs. [H [21 and |3] the two-body spectrum in the trap is plotted for s-wave 
interactions and various values of the background parameters of the resonance, a^^ and 
TgQ, as funciton of a^{B)/alg. We plot both the full solution and the standard case with 
= for comparison. In Fig. [H alg/h = 0.01 which is very small. This means that 
the spectrum is almost equal to the non-interacting case when r^Q is also small. Note 
that the molecular state has energy proportional to — (a^g)"^ so its energy is below the 
range of the figure in the case of r^Q = 0. In Fig. [U^c) and (d) the situation changes and 
a molecular state can be seen. One also clearly see the Zeldovich rearrangement effect 
[22l [23] of the levels in the right of (c) and middle of (d). The connection between the 
Busch model and this effects in the r^Q = case was discussed recently by Farrell et 
al. [23]. However, the a'^{B) value of the rearrangements can now depend on the level 
since we have the term linear in E in Eq. [291 In contrast, for r^Q = the rearrangments 
happens at a'^{B) = 0. 

In Figs. [21 and [3l we exhibit the spectrum for larger values of al^/b, which means 
that the molecular state is now seen even for small r^Q. From these figures it is also 
clear that for large rl^/b the rearrangement happens when a^{B) = a^g. This can 
again be understood from Eq. [23 since x — )■ oo when a^{B) — t- a^^ and which implies 
1 = algf{E)/b. The presence of the linear E term for r^g > gives distortion to 
this simple picture and enriches the rearrangement effect. Notice that in Fig. [3]^a) the 
lowest state shown is in fact not the molecular state but the first excited state. As rfn 



Two-Body Spectra of Harmonically Trapped Atoms 



10 




3 



















f 


(c) 

( 

1 





























(b) 




r 

































(d) 

1 

1 





a={B)/a6., 



Figure 3. (Color online) Same as Fig. [T] for a^g/b — 1 and (a) r^o/b 
rlJb - 1, (c) rlJb = 10, and (d) rljb = 100. 



0.1, (b) 



increases the molecular state is pulled up in energy as seen in Fig. [3]^b). Also, panel (d) 
demonstrates that for very large values of r^Q, the region where the the levels rearrange 
can become very small, yielding almost abrupt jumps in the spectrum. 

To access the spectra above experimentally will require large values of alg/h and 
rlQ/h. For typical traps with h of order /xm, this seems inaccessible. However, a single 
site of an optical lattice could have a much smaller h and has been used before to probe 
the two-body spectrum [6]. Assuming that one could make large reduction in 6, we 
still require Feshbach resonances with large background parameters. Atoms like ^^Na 
or ^''Rb do in fact have known Feshbach resonances that are extremely narrow [8] and 
will give of order /im. Resonances between two different mass atoms also tend to 
be narrow in general so mixtures is an option. However, narrow resonances require the 
ability to tune the magnetic field extremely precisely. The level of tunability required 
here is probably beyond any current experiment but will perhaps be available in next 
generation experimental setups. 



4. Two Dimensions 



The two-dimensional case is similar but contains the peculiarities of 2D scattering 
The free 2D Green's function with out-going boundary condition is 



(30) 



where Hm\x) is the m'th order Hankel function of the first kind. We use the partial 



wave decomposition of the Green's function (see Appendix A ) 

oo 

G%{r,r') = zn ^ J\^\{kr')Hl^^{kr), 



(31) 
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for r > r'. Jm{x) is the Bessel function of order m. The couphng constant, am, has to 
be modified shghtly to fit the 2D geometry. We define 



OLr. 



drJim\ikr)e 



imdr 



(32) 



where 4>{r) = (f){r)W{r) just hke in the 3D case above and 6r is the angle of the 2D vector 
r. For a wave function with angular momentum m, the Fourier transform becomes 

^ CO 

0(fc) = ^ ^"''^'"''"me*"'^'=. (33) 

* m=—oo 

The finite part of the 2D Green's function has to fulfill 



fm{k) 



drdr'(f)*{r)GUr,r')(f){r'), 



(34) 



where G§{r,r') satisfies once again Eq. [13] and the scattering amplitude, fm{k), is 
connected to the scattering phase shift, Sm,{k), through 



f"^W plHcot5^(A;) -iPH- 
This can be solved similarly to the 3D case by assuming that (r > r') 
G'i(r,r') = G|^|(r,r')J|™|(A;r')e*™'''-^"'-', 

which yields 

' f \ \m\ 



GH(r,r') = A{r')U[-Um, \m\ + l^r^b'] 



e 262 



(35) 
(36) 
(37) 



where we define z/^ through E = fvjj{2vm + |^| + !)• Here we have taken a small 
shortcut by introducting the Tricomi hypergeometric function, U{a,b,z), which is the 
convergent solution for z ^ 1. Demanding that (r, r') be regular at the origin yields 
the condition 

2Hr[-i.^] 



lim A{r') 

r'-i>0 



k\" 



(38) 



Proceeding with general m is not attractive since the expressions for the lowest order 
terms are cumbersome. We therefore specialize to specific m values. 

For m = 0, the spectrum should be universal according to the pseudopotential 
approach [3l |35]. The Green's function approach has been discussed in a quasi-2D 
geometry with a tight transverse confinement [36]. First consider the behavior of 
GQ{r,r') at the origin 

'kb' 



Go(r,r') -in + 21n 



^<2-T 



(39) 



for r, r' -> (r > r'), where 'il){x) is the digamma function. The scattering phase shift 
for m = in 2D can be written 



cot 5o(^) = — I 7 + In 



ka 



2D 



+ :^irirk' 



(40) 
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Figure 4. (Color online) Two-body spectrum in 2D for particles interacting in the 
s-wave channel. Solid (black) line is for Tg/fo = 0, dashed (red) for r^/b — 0.5, and 
dotted (blue) has Ve/b — 1.0. 



where is the 2D scattering length and is the effective range. 7 is Euler's constant. 
We have included the effective range term to discuss its effects below. Using Eq. |31l we 
arrive at the eigenvalue equation for the spectrum 

In the case where = 0, this agrees with earlier work [3l |36l [351 EH EH]- In the limit 

)■ 00, the energy approaches the universal expression E = —h^e~'^'^/2ix[a2jjY 

rl = 0) and represents the two-body bound state energy in the absence of the trap 
[40] . This is reasonable since the trap becomes irrelevant for large binding energy and 
small bound state size. The spectrum is shown in Fig. IHfor different values of r*. We 
can see that effective range corrections will alter the energetics of the lowest state quite 
severely within this model. This is very similar to what is found for p-waves with range 
corrections in 3D |20j . 

The Green's function method is particularly transparent when including higher- 
order correction terms in comparison to the pseudopotential approach [3] or, 
equivalently, the Bethe-Peiels boundary condition [HI [38] . A mathematical formulation 
of pseudopotential approaches in any dimension and for any angular momentum 
was recently discussed by Stampfer and Wagner 02j which details the intricate 
problems of even dimensions in comparison to odd dimensions. The expression for 
the pseudopotential beyond lowest order is, however, involved. The Green's function 
approach accomplishes these corrections in a simple manner. 

Whereas the case of m = was universal in the sense that the dependence on 
the details of the potential (through 0(r) dropped out, there are a lot of indications 
that the case \m\ > is not so simple and will depend on short-distance physics. This 
has been pointed out by a number of authors, and leads to the introduction of energy- 
dependent scattering lengths (43], [30l [201 EH] (see [44] for an alternative approach to 
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(41) 
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pseudopotentials that can be applied in a well-defined manner in momentum space in 
both 2D and 3D). We now address this issue within the Green's function approach for 
m = 1. 

The expansion of Gi{r,r') 



IS 



Gi(r,rO 



Tckr 



. 2, 


'kb' 


— z H m 




TT 


_Y_ 



+ 



n{kby 



71 



(42) 



Here we have extracted a factor in front which agrees with Ji{kr) to lowest order in 
kr. We find no other terms that depend explicitly on r. In the m = 1 channel it 
seems natural to define the phase-shift relation that generalizes the m = result as (see 



Appendix B[) 



cot 5i(A;) = — ( 7 + In 

TV 



krjP 



A 



{ka 



2D) 



Ztt 



(43) 



to order fc^ with A a dimensionless constant. The non-universal information about the 
two-body interaction potential is in fact carried by A as pointed out in Ref. [15] and 
discussed further in Appendix B Let us consider the case when A = which occurs if 
there is a bound state at zero energy in the m = 1 potential. In this case we arrive at 
the very simple equation for the eigenspectrum 

"2 1 

+ + l ^ = ln — • (44) 



7 + ^^(-^i) + ^ 



4z/i+4 IU.2D. 

The expression we get is very similar to the m = case, expect for the last term on the 
left-hand side which is the new piece. This extra term will vanish for large energies, but 
will be important around zero energy (where a bound state in free space resides). The 
spectrum is plotted in Fig. [51 We can see that the levels are pushed down compared to 

the m = case, and a state reside at zero energy in the non-interacting -4 )■ limit. 

A bound state appearing below the m = ground state for > is not uncommon 

and occurs also for p- waves in 3D [20]. In the case of A < 0, there is a bound state 
with finite binding energy in the potential. However, we have checked that this only 
gives minor quantitative changes compared to the A = in Fig. [5l Since A = and 
A < yield qualitatively the same spectra, we conclude that the extra term proportional 
to A;~2 makes little difference when the potential has a bound state (possibly at zero 
energy). Note that we are considering the zero-range limit for the potential, so there 
can be at most one bound state. 

The m = 1 spectrum in a 2D trap was discussed in Ref. [35] and a spectrum can 
be found in Fig. 1 of Ref. [3S]. We find good agreement with that result in the case 
when a bound state is present in the two-body potential, corresponding to A < in 
Eq. HHl The level rearrangement takes place at a^^ — )■ as seen in Fig. [5] and there 
is a visible distortion of the manner in which this occurs similar to the examples in 
3D of Fig. [H [21 and [31 These features can also be seen in Ref. [35]. When A > 
in Eq. [131 we find that lowest state seen in Fig. [5l does not diverge to minus infinite 
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Figure 5. (Color online) Two-body spectrum in 2D for p-wave interactions. Dashed 
(red) line is p-wave with r^/b = 0. For comparison, the solid (black) line is s-wave 
with rl/b = 



in binding energy for afj) — >■ 0, but rather behaves similar to the higher lying states. 
This is consistent with the finding in Ref. [SS], although we caution that the scattering 
area used to parameterize the strength in Ref. [35] can have both signs while our a^^y is 
defined to be positive. 

To further explore the dependence of the m = 1 spectrum on the potential 
parameters, we show in Fig. |6] results for a phase-shift of the same form as the hard- 



sphere potential ( [Appendix Bj ) but with positive scattering length so that a bound state 
occurs. This is a somewhat unphysical potential but it helps illustrate the point that 
the spectrum is quite robust under changes in the value of A (for A < 0) since Fig. [5] 
has A = while Fig. [6] has A = —A/it (see Eq. IB.7P . Using a square well instead yields 
almost identical results and we have not plotted this case. In fact, the procedure of 
using a model potential to fix the phase-shift used here is similar to the self-consistent 
energy-dependent pseudopotential methods employed for waves in 3D in Ref. [20] . 
This provides a significant improvement over the energy-independent pseudopotential 
in that case. 

In the case of \m\ > 1 we find that there are terms in G\m\{r, r') that do not vanish 
in the limit of r, r' — t- and that are not represented on the left-hand side of Eq. |34l We 
therefore conclude that the higher partial waves in 2D also yield non-universal spectra. 



5. Conclusions and Outlook 



Using Green's function techniques we have derived the spectrum of two particles in an 
isotropic harmonic trap interacting through a potential that has a range that is much 
smaller than the trap length scale in both three and two dimensional space. In the 
3D case, the spectrum is universal, i.e. independent of the short-range details of the 
interaction, only for the s- and p-wave channels. For rf-wave and beyond there is a 



Two-Body Spectra of Harmonically Trapped Atoms 



15 
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Figure 6. (Color online) 2D p-wave eigenspectrum using a phase-shift similar to the 
hard-sphere potential (solid (black) line). For comparison, the spectrum using the 
generic phase-shift in Eq. |43] is shown by the dashed (red) line. 

dependence on the short- distance physics. We derive a general criterion for when the 
universal expression is valid for higher angular momenta which implies that one must 
be away from any resonances where the scattering length diverges. In the 2D case, only 
the s-wave spectrum is truly universal, but we demonstrated that the p-wave spectrum 
is not very sensitive to the exact details of the interaction as long as it can accomodate 
a two-body bound state. 

The spectra in both 2D and 3D can be interpreted very nicely in terms of Zeldovich 
rearrangement, which occurs when adding a short-range (two-body interaciton) to a 
long-range potential (trap) and tuning through a resonance of the short-range part. 
Including effective range corrections in the 3D s-wave channel shows that one can get 
a very rich set of rearrangment points depending on the background parameters of the 
resonance in the short-range interaction. In the realm of cold gases, this requires use of 
very narrow Feshbach resonances or very tight trapping conditions, both of which are 
currently beyond experimental capabilities but hopefully could be explored in future 
generation experiments. 

It would be interesting to extend the current formalism to polar molecules for which 
external trapping potential are of course also always present in experiments. A number 
of recent works [?6] have explored the bound state structure of such systems but it not 
clear how much influence a harmonic trap or an optical lattice has on these few-body 
states. 
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Appendix A. Free Green's Functions 

The free Green's functions with out-going wave boundary condition used in this paper 
are taken to satisfy the equation 



G%=^-^6{r-r'), (A.l) 



2/i 

where we define the corresponding momentum through E = /i^/c^/2/i. 

Appendix A.l. 3D 
The solution in 3D is 

„ik\r—r\ 

The partial wave expansion that we employ is written in terms of spherical Hankel 
functions of the first kind, h[^\x), and spherical Bessel functions, ji{x). For r > r' we 
have 

G°^(r,r') = AmkJ2Mkr)hf\kr)YUr)Y;^{r'), (A.3) 

Im 

and for r < r' one interchanges the radial variables. 
Appendix A. 2. 2D 

The solution in 2D has been given by many authors, here we follow Adhikari [31] . With 
out-going wave boundary condition, the solution is 

GUr,r') = -^^Hi,'\k\r-r'\), (A.4) 

where Hm\x) is the first Hankel function of order m. The partial wave expansion of 
this Green's function appears to be less accessible and we therefore write it explicitly 
here. Starting from momentum space we define G% = [E — Hq)^^. This implies that 

1 r pip{r-r') 

^""^"""'^ = (2^ y ^'^:^r^' ^^-5) 

where p denotes the wave vector to be integrated over. Using the expansion 

oo 

gifcxcos<A_ zl"lj|„l(A;a;)e^"^'^, (A.6) 

m=— oo 

where Jm{x) is the Bessel function of order m, we can turn this into 

_J^y^ n ^pp '^\m\ipr)J\m\{pr') ^^^^ 

Tlh? In — k"^ 

This integral can be found in standard tables |17] and through analytical contiuation 
we obtain for r > r' 

-^E^H(fcO<|(A:r). (A.8) 

m 
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With our normalization of the Green's function in Eq. lA.H we finally end up with 

GUr,r')=tnJ2JHikr')Hi;;;^ikr). (A.9) 



Appendix B. 2D Scattering 

Scattering in 2D is complicated by the appearance of logarithmic terms in the typical 
wave function in the asymptotic region of large distance, which is a Neumann function, 
Ym{x). Note that we are only interesting in short-range potentials (vanishing for 
distances r > r^) for which the asymptotic solution is the free one. We can therefore 
write the angular momentum m scattering wave function, for r > tq in the form 

^„(r) = A [cot 6m{k)Jm{kr) - Ym{kr)] (B.l) 

where E = h'^k'^/2^ and is a normalization constant. The phase-shift Sm{k) can be 
calculated from 



cot 6mik) 



(B.2) 



where x = krQ and prime denotes derivative with respect to x. The logarithmic 
derivative is 

1 c/^^fr)" 



7« 



\1/ 



dr 



(B.3) 



J r=ro 



We would like to discuss this in terms of an appropiately defined scattering length. 
Here we follow the intuitively clear defintion |18] 



"^oir) -> In 
^^fr) r' 



r 

ao 



and 



1 - 



n \ 2r7i 
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(B.4) 
(B.5) 



for m > 0. 

The m = case was studied by Verhaar et al. [37] who found the expression 



cot 5o(^) = — I 7 + In 

TT 



In 



(B.6) 



to second order in k. We expect a similar expression for the m = 1 phase-shift, but 
very little can be found on this in the literature. To check this, we compute the exact 
expression for the hard-sphere and for the attractive square well potential. The hard- 
sphere yields 

4 3 2 

cot h{k) = — - - — + - 

7T{kro) 271 7T 

and the square well gives 



7 + ln 



kro 



(B.7) 



cot 6i{k) 
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in the limit 2firQVQ/h? — )■ where the depth is — Vq. By using the definitions of the 
scattering length, a^, above, we find af = r^ for the hard-sphere and = y for 
the square well. The choice of sign for ai is not given directly by these relations. In 
the regime where there is a bound states in the potential, the definition above implies 
Qm > 0. Since we know that the hard-sphere potential should not hold a bound state, 
a suitable choice is ai = — rg. Similarly for the square well in the limit 2/iroVo//i^ — )■ 
where the centrifugal barrier hinders the formation of a bound state |19]. We will be 
interested in the regime ai > only, as it is nicely comparable to the m = case. To 
approach the zero-range limit with a bound state always present for the square well, 
one needs to take the limit of tq — )■ and —Vq — ?■ oo in a manner that keeps 2firQVQ/ff 
at or above the critical limit for the appearance of a bound state. 

From the discussion above, we get the suggestsive expression for the phase-shift 



A 2 
cot 6i{k) = - — - + 5 + - ( 7 + In 



kai 



(B.9) 



where A and B are potential-dependent low-energy constants. We find that the structure 

of cot5i(fc) is very similar to the m = case, except that there is an added k~^ oc 

term. The leading divergence in the corresponding scattering amplitude is therefore a 

pole rather than the logarithm as for m = 0. This was pointed out in Ref. [15] • In the 

main text, we have shown that in the presence of a harmonic trap, it does not make 

much difference whether the leading term is included or not when calculating the 

two-body spectrum in the presence of a bound state. This can be seen from comparison 

of Fig. [5] (with A = 0) and Fig. |6] (with A < 0). However, in the case with A > there 

is no bound state as discussed in the main text. 
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